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ABSTRACT 

We present results from a WIYN 3.5-m telescope imaging study of the globular 
cluster (GC) systems of the edge-on spiral galaxies NGC 891 and NGC 4013. We 
used the lO'xlO' Minimosaic Imager to observe the galaxies in BVR filters to 
projected radii of ~20 kpc from the galaxy centers. We combined the WIYN 
data with archival and published data from the WFPC2 and Advanced Camera 
for Surveys on the Hubble Space Telescope to assess the contamination level of the 
WIYN GC candidate sample and to follow the GC systems further in toward the 
galaxies' centers. We constructed radial distributions for the GC systems using 
both the WIYN and HST data. The GC systems of NGC 891 and NGC 4013 
extend to 9±3 kpc and 14±5 kpc, respectively, before falling off to undetectable 
levels in our images. We use the radial distributions to calculate global values for 
the total number (Ngc) and specific frequencies {Sn and T) of GCs. NGC 4013 
has Ngc = 140±20, Sn = 1.0±0.2 and T = 1.9±0.5; our Ngc value is -40% 
smaller than a previous determination from the literature. The HST data were 
especially useful for NGC 891, because the GC system is concentrated toward 
the plane of the galaxy and was only weakly detected in our WIYN images. 
Although NGC 891 is thought to resemble the Milky Way in terms of its overall 
properties, it has only half as many GCs, with Nqc = 70±20, Sn = 0.3±0.1 
and T = 0.6±0.3. We also calculate the galaxy-mass-normalized number of blue 
(metal-poor) GCs in NGC 891 and NGC 4013 and find that they fall along a 
general trend of increasing specific frequency of blue GCs with increasing galaxy 
mass. Given currently available resources, the optimal method for studying the 
global properties of extragalactic GC systems is to combine HST data with wide- 
field, ground-based imaging with good resolution. The results here demonstrate 
the advantage gained by using both methods when possible. 

Subject headings: galaxies: formation — galaxies: individual (NGC 4013, NGC 891) 
— galaxies: spiral — galaxies: star clusters: general 
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1. Introduction 

Globular clusters (GCs) are among the most useful stellar populations to study in 
galaxies. These compact star clusters each contain ~10^ — 10^ of stars that appear to have 
formed over a relatively short period of time and under similar conditions (e.g., Ashman & 
Zepf 1992; Ashman & Zepf 1998; West et al. 2004; Brodie & Strader 2006; Griffen et al. 2009; 
Gnedin 2010). GCs are luminous (Milky Way GCs typically have My ~ —10 to —5; Harris 
1996), often many Gyr old (e.g., Chaboyer et al. 1998, Brodie & Strader 2006), appear to be 
associated with major star formation episodes in galaxies (Brodie & Strader 2006, Bastian 
2008), and may be forming in nearby galaxy mergers and/or starbursts (e.g., Whitmore & 
Schweizer 1995, de Grijs et al. 2001). GCs are detected in virtually all giant galaxies that 
have been imaged to sufficient depth to cover a significant fraction of the GC luminosity 
function (GCLF; Harris 1991, Brodie & Strader 2006). Together, these properties mean 
that GCs can be used as a stellar record of the major events associated with the origin and 
evolution of the galaxies that host them. The numbers, spatial distributions, luminosities, 
metallicities, and kinematics of the GCs provide fundamental information about the star 



formation and assembly histories of giant galaxies ( iBrodie fc Strader 1120061 ). 



Motivated primarily by these ideas, we have been carrying out a wide-field CCD imaging 
survey of the GC systems of giant spiral, eUiptical, and SO galaxies in the local universe 
(distances to ~20 Mpc). The spiral galaxies we target are all edge-on or nearly edge-on, 
to make it easier to identify bona fide GCs and less likely that we will mistake knots of 
star formation or open clusters in the galaxy disk for GCs. We use large-format and mosaic 
CCD detectors to image the target galaxies in multiple broadband filters and select as GC 
candidates the point sources around the galaxies with the expected magnitudes and colors. 
We use the GC candidate lists and galaxy images to estimate the global properties of each 
galaxy's GC system — namely, the total number and specific frequency of GCs as well as 
the spatial and color distribution of the overall system. Since GC systems of nearby galaxies 
are extended (often covering tens of arc minutes from the galaxy center) , wide-field imaging 
is necessary to cover enough of the GC system to allow for direct, accurate determinations 
of these quantities. 

The general goal of the survey is to increase the number of galaxies for which accurate 
measurements of the global GC system properties exist, so that these measurements can in 
turn be used to test scenarios for the formation and evolution of giant galaxies. The details 
of the design and methods of the survey, as well as the results from the first ten galaxies 
studied, are given in Rhode & Zepf (2001, 2003, 2004; hereafter RZOl, RZ03, RZ04) and 
in Rhode et al. (2005, 2007; hereafter R05, R07). Another motivation for the survey is to 
use the GC candidate lists as target lists for follow-up spectroscopy, in order to measure 
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the radial velocities of the GCs around each galaxy. The kinematics of the GCs provide 
additional insight into the origin of each host galaxy and the radial velocities can also be 
used to quantify the galaxy's mass distribution (e.g., Bergond et al. 2006, Bridges et al. 
2007). 

In this paper, we present results for two more galaxies imaged as part of the wide- 
field survey: NGC 4013 and NGC 891, both edge-on Sb spiral galaxies. Some fundamental 
properties of the galaxies are summarized in Tabled! To analyze these galaxies' GC systems, 
we combined wide-field mosaic CCD data obtained with the WIYN 3.5-m telescop^ with 
previously published results from Hubble Space Telescope (HST) imaging studies of the inner 
parts of the galaxies' GC systems. Ground-based imaging and HST imaging each have 
advantages and disadvantages as methods for studying extragalactic GC systems. Ground- 
based imaging with large-format CCDs often enables one to image all, or nearly all, of the 
galaxy's GC system in a single pointing. However, limited resolution makes it difficult to 
distinguish GCs from distant background galaxies and to detect GCs close in to the center 
of the parent galaxy, especially in the case of a spiral galaxy with a dusty disk and/or 
active regions of star formation. HST's superb resolution allows one to detect GCs further 
in toward the centers of their host galaxies, to distinguish GCs from compact background 
galaxies, and in some cases to study the sizes and structural parameters of the individual 
GCs (e.g., Spitler et al. 2006, Forbes et al. 2010, Harris et al. 2010). On the other hand, the 
limited fields-of-view of the HST imagers means that only a fraction of the GC system is 
typically observed (see discussions in Ashman & Zepf 1998 and Brodie & Strader 2006), and 
properties such as the overall spatial distribution and specific frequency of the GC system 
can be undetermined or highly uncertain without supplemental data or assumptions (e.g., 
Goudfrooij et al. 2003, Chandar, Whitmore, & Lee 2004, Peng et al. 2008, Cantiello, Brocato, 
& Blakeslee 2009). 

The results presented here for NGC 4013 and NGC 891 demonstrate the efficacy of 
combining wide-field ground-based data with data from HST to study GC systems of giant 
galaxies in the local universe. The following section describes the WIYN observations and 
initial image reduction steps. Section [3] presents the techniques used to detect and select GC 
candidates in the WIYN images. Section H] describes additional analysis steps we performed 
to assess the completeness and contamination level of the WIYN data and ends with a 
description of how we combined the published HST results for these galaxies with the WIYN 
data. The overall results for the galaxies' GC systems are described in Section [51 The paper 
ends with a summary of our main conclusions in Section ^ 



^The WIYN Observatory is a joint facility of the University of Wisconsin, Indiana University, Yale 
University, and the National Optical Astronomy Observatory. 
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2. Observations and Image Reductions 

NGC 891 and NGC 4013 were imaged in ttiree broadband filters (BVR) witli tlie Mini- 
mosaic camera on the 3.5-m WIYN telescope at Kitt Peak National Observatory. The WIYN 
Minimosaic is constructed of two 2048x4096-pixel CCDs with 0.14" pixels and the field-of- 
view is 9.6' on each side. NGC 891 is relatively nearby (see Tabled]) and thus was positioned 
on the edge of the detector to maximize the radial coverage of its GC system; NGC 4013 is 
more distant and was placed in the center of the field. Figures [1] and [2] show the locations 
of the WIYN fields, superposed on gray-scale images from the Digitized Sky SurvejQ. The 
WIYN observations analyzed for this paper were acquired in 2001 January and 2009 March; 
Table [2] lists the dates, number of images, and exposure times used for each galaxy pointing. 

The sky conditions were photometric during the second night of the 2001 January ob- 
serving run. On that night, images of standard stars from Landolt (1992) were taken, along 
with one or more exposures of NGC 891 and NGC 4013, in all three filters. A set of pho- 
tometric calibration coefficients (color terms and zero point constants) was derived using 
the Landolt star observations taken that night. The color coefficient in the V magnitude 
equation was 0.02±0.01. The color coefficients in the B — V and V — R color equations 
were 1.02±0.01 and 1.05±0.03, respectively. The formal errors on the zero points in the V, 
B — V and V — R calibration equations ranged from 0.004 to 0.006, confirming that the sky 
conditions that night were indeed photometric. 

Images of the target galaxies taken on other, non-photometric nights during the 2001 
January run were calibrated by scaling them to a selected image taken on the photometric 
night and then combining them with that image (see below for details of the combining 
process). For NGC 4013, a set of V images was obtained in 2009 March; those images 
were postcalibrated by calculating the zero-point offset in V magnitude between the final, 
combined 2009 March V image and a single V image taken in 2001 January. 

Initial reductions (overscan and bias level subtraction, flat-field division) were performed 
on the Minimosaic images using standard reduction tasks in the IRAF package MSCRED. 
The tasks msczero, msccmatch, and msc image were then used to transform the multi- 
extension FITS images into single-extension images. All of the images for a given galaxy 
were aligned to match a single pointing. A constant sky background level was computed and 
then subtracted from each individual galaxy exposure. Images of the same galaxy target. 



^The Digitized Sky Surveys were produced at the Space Telescope Science Institute under U.S. Govern- 
ment grant NAG W-2166.The images of these surveys are based on photographic data obtained using the 
Oschin Schmidt Telescope on Palomar Mountain and the UK Schmidt Telescope. The plates were processed 
into compressed digital form with the permission of these institutions. 
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taken in the same filter, were scaled to a common flux level and combined using the IRAF 
task imcombine and the ccdclip pixel rejection algorithm. The constant sky background 
level measured in the image used as a reference in the scaling step was added to the final 
combined image. The end result was one deep, stacked image per filter for each galaxy 
target. The mean full-width at half-maximum of the point-spread function (FWHM PSF) in 
the final combined images ranges from 1.0"-1.2" for NGC 891 and 0.7"-0.8" for NGC 4013. 

3. Detection and Analysis of the Globuleir Cluster System 

3.1. Source Detection 

GCs should appear in our ground-based WIYN Minimosaic images as a slight over- 
density of point sources surrounding the target galaxy. To detect and then analyze the GC 
systems of NGC 891 and NGC 4013, we executed the same series of steps used in the previous 
studies that are part of our wide-field survey (RZOl, RZ03, RZ04, R07). We first smoothed 
the combined images using a circular median filter with a diameter equal to ~7.25— 7.5 times 
the mean FWHM PSF of the image. The smoothed images were subtracted from the original 
versions to remove the diffuse light from the target galaxy and make it easier to locate the 
point-source GCs. The constant sky background level present in the original images was 
added to the galaxy-subtracted images; all subsequent analysis steps were performed on the 
galaxy-subtracted images. The IRAF task daof ind was used to detect sources above a given 
signal-to-noise level in each frame. Areas of the images likely to produce spurious source 
detections - i.e., near saturated stars and in the dust-obscured spiral disks of the target 
galaxies - were excluded from the detection process. A final fist of sources detected in all 
three filters was produced; there were 774 objects in NGC 891 and 998 in NGC 4013. 

3.2. Removing Extended Sources 

Objects that appear extended in our ~l"-resolution images are likely to be background 
galaxies rather than GCs. To remove these objects, we measured the FWHM of each source 
and plotted it versus the source's instrumental magnitude. Bright point sources form a tight 
sequence around the mean FWHM value for each image, with increasing scatter around the 
mean as the instrumental magnitude increases. We wrote software to select objects based on 
this pattern; the range of acceptable FWHM values for point sources increases with source 
magnitude. If a source failed the FWHM criteria for a point source of a given magnitude in 
any one filter, it was removed from the list of potential GC candidates. After this step, the 
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number of sources remaining was 490 for NGC 891 and 488 for NGC 4013. 



3.3. Aperture Photometry 

Aperture photometry was performed on the remaining point sources in the BVR images 
of NGC 891 and NGC 4013, using an aperture of radius equal to the mean FWHM PSF for 
the specific image. An aperture correction was derived for each image by measuring the mean 
difference in the total magnitude and the magnitude within a one-FWHM-radius aperture. 
The aperture corrections ranged from —0.229 to —0.308 with errors of 0.001 to 0.005. Final, 
calibrated total magnitudes were calculated for the point sources in the images by applying 
the appropriate aperture correction and photometric calibration coefficients, as well as a 
Galactic extinction correction derived from the reddening maps of Schlegel, Finkbeiner, & 
Davis (1998). The Galactic extinction corrections for NGC 891 were Ab= 0.277, Ay=0.209, 
and Ar=0.168, and for NGC 4013 were Ab= 0.073, Ay=0.055, and A/{=0.045. 



3.4. Final Selection of GC Candidates 

The last step in the GC candidate selection process was to choose objects with BVR 
magnitudes and colors that match those expected for GCs at the distance of the host galaxy. 
We executed this step using the same basic method described in the previous survey papers 
(RZOl, RZ03, RZ04, and R07). GCs were assumed to have My^ —11 to —4, based on the 



GCLF of the Milky Way and other giant galaxies (j Ashman fc Zepf Ill998l ). Objects within 
this magnitude range (given the galaxy distances in Table [1]) were checked to determine 
whether their BVR colors were like those expected for GCs with metallicities between [Fe/H] 
of —2.5 and 0.0. Specifically, sources with B — V and V — R values (and corresponding 
photometric errors) that placed them within 2— 3-cr of the derived relation for Galactic GCs 
with the specified metallicity range were accepted as GC candidates. 

As explained in R07, the spiral galaxies included in our wide-field survey typically have 
so few GC candidates (a few dozen, compared to hundreds or thousands for the elliptical 
galaxies) that we carefully examine each of the GC candidates as well as those that narrowly 
miss meeting the above criteria. We then adjust the selection criteria accordingly, so as 
not to exclude objects that are likely to be GCs or mistakenly include objects that are 
probably contaminating foreground stars or background galaxies. We consider things like 
where objects lie in the BVR color-color plane, where they are located relative to the galaxy, 
and how the objects appear in archival i/S" T images (see Section [4. 2. 2p . 



- 7- 



For NGC 4013, we executed the usual selection steps as described above, accepting 
sources with V magnitudes >19.9 and BVR colors within 3-cr of the B — V vs. V — R 
relation for Galactic GCs. We also decided to impose a faint magnitude limit oiV < 24.4, or 
~0.7-(j past the GCLF peak, in order to minimize the contamination from faint background 
galaxies that were not removed in the extended source cut. The final set of GC candidates 
in NGC 4013 includes 83 objects. 

The GC candidate selection for NGC 891 was much more involved. This galaxy is at 
relatively low Galactic latitude (—17 degrees vs. the more typical range for our targets, +40 
to +70 degrees) and therefore has a higher surface density of Galactic stars than typical 
galaxy images from the survey. Imposing the usual luminosity threshold of My = —11 
yielded a significant population of bright sources in the GC candidate list when, based on 
the shape of the Milky Way GCLF and that of other giant galaxies, only a few real GCs 
would be expected at those magnitudes. To account for this and to limit contamination 
from background galaxies, we imposed more stringent V magnitude criteria on both ends of 
the luminosity function, choosing GC candidates in the range 20.7 < V <23.6 (extending 
from ~1.3-(T brighter than the GCLF peak to l-cr fainter than the peak). We also used a 
2-(j criterion for the BVR color selection because analysis of HST images indicated that 
constraining the color cut would help eliminate contaminants from the sample, while still 
allowing real GCs to be selected. Finally, we explicitly excluded 13 objects that had count 
ratios in HST images that indicated that they were background galaxies or foreground stars. 
(The HST data analysis is discussed in Section r4.2.2p . 

Lastly, we added a few sources to the final hst for NGC 891. We included four sources 
that were located very near the galaxy disk and had colors that put them just outside the 
color selection box, in the direction of the reddening vector. We added two more sources that 
had been excluded by the more stringent V magnitude criteria, but nevertheless appeared to 
be real GCs based on the i/ST study we used to supplement the WIYN data (see Section 1^^ . 
The final set of GC candidates for NGC 891 includes 43 objects. 

Figures [3] and H] show the results of the color selection. The point sources in each 
galaxy field that appeared in all three filters are marked with open squares and the final 
GC candidates are marked with filled circles. The 2- or 3-cr-wide color selection boxes are 
also marked. Note that due to the V magnitude criteria, not all objects that appear within 
the color selection boxes are GC candidates. For illustrative purposes, tracks are plotted in 
Figures [3] and m that show where galaxies of different morphological types (E/SO through 
Irregular) at different redshifts would lie in the color-color plane. RZOl details how the tracks 
were calculated. The distribution of sources looks markedly different in the NGC 4013 and 
NGC 891 fields and the galaxy tracks help explain why. The NGC 4013 images are both 
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deeper (see Section 14.11) and located at higher Galactic latitude (+70 degrees) than the 
NGC 891 images. Consequently the figure is dominated by a large number of objects located 
in the lower left part of the BVR color-color plane for NGC 4013; these are most likely faint, 
blue, low- to mid-z background galaxies that are unresolved in our ground-based images. By 
contrast the NGC 891 field has relatively few objects in the lower left corner of the BVR 
plane, and instead has a very well-populated, narrow sequence of objects running diagonally 
through the plane; these are likely Galactic foreground stars. 

The color-magnitude diagrams (CMDs) of the GC candidates for each galaxy are shown 
in Figure |5l To make use of information from all three filters, the V magnitudes vs. the 
B — R colors of the sources are plotted. The two very red sources with B — R ~1.9 in 
the NGC 891 CMD are very close to that galaxy's disk, which may help explain their red 
colors; the brighter of those two sources, with ~ 19, was added to our sample based on its 
appearance as a legitimate GC in an HST ACS study from the literature (see Section 1^^ . 

4. Additional Analysis Steps 

4.1. Completeness Testing and Detection Limits 

We quantified the point-source detection limits of the WIYN Minimosaic images by 
performing a series of completeness tests. In each test, 200 artificial sources, each with a 
point-spread-function that matched the best-fit PSF for the image and a total magnitude 
within 0.2-magnitude of a specified value, were added to a given Minimosaic image. The 
same detection steps used on the original images were then executed and the fraction of ar- 
tificial sources that was recovered was recorded. The process was repeated at 0.2-magnitude 
intervals over a range of 5—6 magnitudes per filter. This yielded a measurement of com- 
pleteness vs. magnitude for each filter and each galaxy field. The NGC 891 images are 50% 
complete at B =25.2, V =24.2, and R =24.1. The images of NGC 4013 were taken under 
better seeing conditions and accordingly are somewhat deeper, with 50% detection limits at 
B =25.6, V =25.2, and R =24.8. 

4.2. Contamination 

A serious concern for ground-based imaging studies of extragalactic GC systems is con- 
tamination from point-source stars and compact galaxies masquerading as GCs in the can- 
didate lists. We went through a series of analysis steps to try to both quantify and correct 
for the contamination from stars and galaxies in the GC candidate lists. 
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4-2.1. Model Prediction for Stellar Contamination 

We used a Galactic structure model (Mendez & van Altena 1996, Mendez et al. 2000) 
to help us evaluate the level of stellar contamination that exists in the GC candidate lists. 
Quantities like the position of the Sun within the Galaxy and the fraction of stars in the 
Galaxy's structural components (disk, thick disk, and halo) are input to the model; it then 
predicts the surface density of stars within a given magnitude and color range and direction 
on the sky. Our experiments with the model indicate that the predicted stellar surface 
density does not strongly depend on the parameters like the Galaxy composition and the 
solar distance and height above the Galactic disk. The model predicts that the surface 
density of Galactic stars with V magnitudes and B — V colors in the same range as the 
samples of GC candidates is 0.06 arcmin"^ for NGC 891 and 0.07 arcmin"^ for NGC 4013. 
Note that the value for the NGC 891 field is fairly small, even though Figure H] suggests a 
large number of stars in the field; this reflects the fact that we specifically reduced the stellar 
contamination in the GC candidate sample for that galaxy by limiting the sample to 20.7 
< V < 23.6. (If we instead impose the usual V magnitude threshold of My = —11, the 
V range becomes 18.6 < V < 23.6 and the stellar contamination predicted by the model is 
nearly four times larger, at 0.23 arcmin"^.) 

4.2.2. Examining the WIYN GC Candidates in Archival HST Images 

Next we investigated the level of contamination in the GC lists due to background 
galaxies. HST can resolve many compact background galaxies that appear as point sources 
in ground-based data, so we analyzed archival HST images of NGC 4013 and NGC 891 to 
determine whether any of the WIYN GC candidates were actually galaxies. 

For NGC 4013, we retrieved Wide Field Planetary Camera 2 (WFPC2) images from the 
HST data archive). We analyzed two pointings in F450W located ~0.5' from the galaxy 
center [HST program GO. 8242, PI: Savage) and one pointing in F814W located ~0.7' from 
the galaxy center [HST program GO. 6685, PI: Huizinga). We requested that "On-the-fly" 
processing be applied to the images before retrieval and we used the STSDAS task gcombine 
to stack multiple pointings taken in the same filter. Twenty-four of the WIYN GC candidates 
appear in one or more of the WFPC2 frames. Photometry was performed on these objects 
using apertures 0.5- and 3-pixels in radius and the ratio of the flux within the two apertures 



■^Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the data 
archive at the Space Telescope Science Institute. STScI is operated by AURA, under NASA contract NAS 
5-26555. 



- 10 - 



was computed. The ratio of counts in the larger aperture to smaller aperture should be 
significantly larger for extended background galaxies than it is for compact GC candidates 
(Kundu et al. 1999). This test indicated that none of the 24 WIYN GC candidates was a 
background galaxy. The 24 candidates were inspected visually to confirm the results of the 
count ratio test. Note that this does not necessarily mean that there is zero contamination 
from galaxies in the WIYN sample, just that the specific subset of candidates (~30% of the 
total) in the WFPC2 images does not appear to include any resolved galaxies. 

For NGC 891, we downloaded WFPC2 images from the HST archive and Advanced 
Camera for Surveys (ACS) images from the Hubble Legacy Archive (HLA)E|. The WFPC2 
images we used were: one pointing in F814W at 0.2' from the galaxy center {HST pro- 
gram GO. 9042, PI: Smartt), a pointing in F450W at r~3.5' {HST program GO. 8805, PI: 
Casertano), and a pointing at r~8' observed in F450W and F606W {HST program 9676, 
PI: Rhoads). These WFPC2 data were processed and analyzed in the same manner as de- 
scribed above. Thirteen of the WIYN GC candidates were located in one or more of the 
WFPC2 frames, and photometry with 0.5-pixel and 3-pixel apertures indicated that one 
candidate was a galaxy. The ACS images we analyzed consisted of three pointings in the 
halo of NGC 891 that overlap the region covered by the WIYN pointing. All three pointings 
were taken in both the F606W and F814W filters and all were from HST program GO. 9414 
(PI: de Grijs). The pointings were located 2.1', 3.7', and 6.8' from the galaxy center. Im- 
ages from the HLA are processed using the standard HST pipeline and combined using the 
Multidrizzle software. We measured the fiuxes in two different apertures for 21 of the 
WIYN GC candidates in the ACS images and found that three of these were galaxies and 
ten were foreground stars. As a result of this analysis, we decided to constrain the GC can- 
didate selection applied to the WIYN data by using a 2-a color selection and by explicitly 
excluding some of the background galaxies or stars that remained even after the 2-a criterion 
was applied. 

Our ability to draw general conclusions about the contamination level in the WIYN 
data was hampered by the small numbers of objects we could identify and measure in the 
HST images. Therefore, rather than using the statistics from the HST analysis (e.g., one 
galaxy out of 13 GC candidates in the WFPC2 images of NGC 891) to apply some sort 
of global contamination correction, we opted to simply remove the extended objects from 
the WIYN GC candidate list and use another method to calculate an overall contamination 



^Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the Hubble 
Legacy Archive, which is a coUaboration between the Space Telescope Science Institute (STScI/NASA), the 
Space Telescope European Coordinating Facility (ST-ECF/ESA) and the Canadian Astronomy Data Centre 
(CADC/NRC/CSA). 
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correction. 



4-2.3. Contamination Correction Based on the Asymptotic Behavior of the Radial Profile 

We constructed an initial radial profile for each galaxy's GC system by assigning the 
GC candidates to I'-wide concentric annuli centered on the galaxy nucleus. We computed 
the effective area of each annulus (the total area minus the regions of the images that were 
excluded from the daof ind search described in Section [3]) and used that to calculate the 
surface density of GCs for the annulus. GC system radial distributions that are created 
from wide-field images in this manner typically exhibit a peak surface density of GCs near 
the center of the galaxy and then fall off monotonically with increasing radius. At some 
distance from the galaxy center, they then flatten out to a fairly constant surface density, as 
long as the images are large enough to have covered the full radial range of the GC system 
(e.g., Harris et al. 1985, Harris 1986, RZOl, RZ04, R07). The constant surface density in 
the outermost annuli can be assumed to be due to contaminating objects, i.e., point sources 
with magnitudes and colors like GCs that cannot be easily differentiated from real GCs in 
a photometric study. Thus the surface density in the outer bins can be subtracted from the 
overall radial distribution to correct it for contamination. 

There were factors that complicated applying this type of asymptotic correction for both 
NGC 4013 and NGC 891. NGC 4013 has recently been shown to have a tidal stellar stream 



( jMartinez-Delgado et al. 1 120091 ) . which appears as a low-surface-b rightness loop extending 



6^ (^ 26 kpc) out from and around the northeast side of the galaxy. iMartinez-Delgado et al. 



( I2OO9I ) analyzed the location of the stream relative to orbital model predictions and concluded 
that it was likely due to a dwarf satellite of mass ~6 x 10^ solar masses that is merging with 
NGC 4013. The highest-surface-brightness portion of the loop appears in our Minimosaic 
images on the northeast side of the disk. There may be a slightly enhanced GC population 
on that side of the galaxy: when one compares the location of the stream relative to the 
GC candidates we detect, it does appear that the candidates might be preferentially located 
within the stream, although the effect is very subtle and may not be real (see Section I5l2|) . 
It does seem possible that the dwarf galaxy responsible for NGC 4013's stellar stream could 
have GCs associated with it, based on our current picture of the best-known accreting dwarf 
galaxy, the Sagittarius dwarf spheroidal galaxy (Sgr dSph). Sgr dSph is being tidally dis- 
rupted by its interaction with the Milky Way and studies of its orbit and stellar populations 
indicate that its mass is ~10^ M0 and it may have deposited as many as eight GCs into the 
Milky Way GC system (Ibata, Gilmore, & Irwin 1994; Johnston et al. 1999; Layden & Saraje- 
dini 2000; Cohen 2004; Carraro et al. 2007; Belokurov et al. 2007; Carraro 2009). In order to 
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avoid over- correcting for contamination (and thus under-counting GCs) in NGC 4013's GC 
system, we therefore excluded 10 GC candidates that might be associated with the stream 
from the asymptotic contamination correction. With these objects excluded, the radial pro- 
file of GC candidates is fairly constant in the outer two (of five) annuli and the mean surface 
density in these annuli is 0.38±0.09 arcmin"^. (Leaving in the 10 GC candidates leads to 
a mean surface density in the outer two annuli of 0.62±0.12 arcmin"^. Note also that the 
10 GC candidates that might coincide with the stream were removed only for the purpose 
of calculating an asymptotic contamination correction and were included in all subsequent 
analysis steps.) 

NGC 891 posed a problem for other reasons. The initial radial profile created by assign- 
ing the GC candidates to 1-' bins did not have the usual shape for a GC system, but instead 
was relatively flat over the entire radial extent of the WIYN data. The surface density of GC 
candidates fluctuated between 0.41—0.96 arcmin"^ in the inner five annuli (at r~l'— 5' from 
the galaxy center), before settling down to values between 0.27—0.40 arcmin— 2 in the last 
three annuli (at r~6'— 8'). We assumed that the last three annuli provided a reasonable es- 
timate of the surface density of contaminating objects, so we took the asymptotic correction 
to be the mean of these values: 0.33±0.04 arcmin"^. 

The asymptotic values derived from the initial radial profiles were adopted as the final 
contamination corrections for NGC 4013 and NGC 891. We note that the stellar contam- 
ination levels estimated from the model (Section 14.2. ip were in both cases lower than the 
asymptotic contamination corrections, which is consistent with the idea that additional con- 
tamination comes from unresolved galaxies. A contamination fraction was calculated for 
each annulus in the radial profile by multiplying the surface density of contaminants by the 
effective area of the annulus and then dividing that number by the total number of GC 
candidates in the annulus. This radially-dependent contamination correction was later used 
to correct the GCLF (Section 14. 3p and the final radial profile (Section 15.10 . 

4.3. Determining the GCLF Coverage of the WIYN Data 

Observed GC luminosity functions were calculated for NGC 891 and NGC 4013 by 
assigning the V magnitudes of the final GC candidates to bins 0.3 magnitude wide. A 
contamination correction was applied to the data by multiplying the number of GCs in a 
particular radial bin by the contamination fraction at that radius, as determined from the 
asymptotic correction. The numbers in each bin were then divided by the total completeness 
correction for that magnitude (the total completeness is calculated by convolving the indi- 
vidual completenesses in the BVR filters, as described in detail in RZOl). The corrected GC 
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luminosity function for each galaxy was then fitted with a Gaussi an function with a pea k 
absolute magnitude like that of the Milky Way GCLF, at My =-7.3 jAshman fc Zepf lllQQSf ). 
corresponding to an apparent magnitude of V= 22.3 and 23.6 for NGC 891 and NGC 4013, 
respectively. We varied the dispersion of the Gaussian between 1.2 and 1.4 mag and found 
that the mean fraction of the theoretical GCLF covered by the observed luminosity function 
was 0.69±0.01 for NGC 891 and 0.656±0.002 for NGC 4013. Changing the bin size of the 
LF data, or excluding one bin from the fit, changed the mean fractional coverage by ~2— 6%; 
we included this uncertainty in the error calculation in Section 15.31 



4.4. Combining HST Results from the Literature with the WIYN Results 

The GC systems of the two galaxies we observed with WIYN have been studied using 
HST data: Goudfrooij et al. (2003; hereafter G03) analyzed NGC 4013 and Harris et al. 
(2009; hereafter H09) analyzed NGC 891. The HST studies are complementary to our 
analysis: they have smaller radial coverage but higher resolution than the ground-based 
data and allow a search for GCs further into the galaxy disks. We used the results from 
these studies in combination with the WIYN data to produce a radial distribution for each 
galaxy's GC system. We include the HST results in all subsequent analysis steps. 

G03 analyzed HST WFPC2 images of NGC 4013 and six other spiral galaxies. The 
images of NGC 4013 were taken in the F555W and F814W filters and cover the central 
~1.5— 2' around the galaxy. The WFPC2 field-of-view is plotted in Figure [21 along with the 
location of the WIYN pointing. G03 used steps similar to those we executed on the WIYN 
images to detect and select point-source GC candidates by their V magnitudes and V — I 
colors. They published a list of the 50 brightest GC candidates in NGC 4013; these have V 
= 21.8 to 24.4. The V^=24.4 faint limit matches the limit we imposed on the WIYN data. 
We compared each object in the G03 list to the WIYN GC candidate list, to examine in 
detail the differences in GC selection in the two data sets. Twenty-five of the 50 objects in 
the G03 list are located in the galaxy disk, in areas of the WIYN images that were masked 
out because we could not reliably detect GCs there. Another 20 of the 50 GC candidates 
appear in both the WIYN list and the G03 list. Of the five remaining from the G03 list, 
three barely appear in our WIYN images and fall below our source detection limits. Two 
others appear in the WIYN images but were eliminated as GC candidates because they had 
large FWHM values or BVR colors outside the expected region. We retained all 50 GC 
candidates in the G03 list to use in the subsequent analysis steps. 

We assumed that the contamination level in the HST candidate list was negligible; 
G03 used detection software that characterized the sharpness and roundness of a source to 
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remove extended objects, and for magnitudes brighter than V^24 the contamination from 
unresolved background galaxies in these types of HST images should be very low (Kundu et 
al. 1999). To estimate the GCLF coverage of the HST data, we plotted the completeness- 
corrected luminosity function of the GC candidates and fitted it with Gaussians of varying 
dispersions, in the same manner as for the WIYN sample. The mean fractional coverage 
(averaged for three different Gaussian dispersions) is 0.75±0.01. The completeness curves 
were taken from G03; the completeness correction was minimal, because the HST data are 
~90-100% complete to V=24A. 

H09 used deep HST ACS data to identify GC candidates in three fields (which they 
designated HI, H2, H3) covering part of the disk and inner halo of NGC 891. The locations 
of the H09 ACS fields, relative to the WIYN pointing, are shown in Figure [TJ Each field 
was imaged in the F606W and F814W filters and H09 selected a sample of possible GC 
candidates via visual inspection. They refined the sample based on the size, shape, ellipticity, 
magnitude, and color of each candidate; their final published list of objects consists of 16 
"Best" GC candidates and 27 more that may be GCs or contaminants. To simplify the 
process of combining the HST and WIYN results, we opted to include in our analysis only 
the sources in the HI and H2 fields from H09. The H3 field covered an area beyond the 
boundaries of the WIYN images and H09 identified only five candidates in the field: two 
"best" GC candidates and three they tentatively label as galaxies; we did nothing more with 
the H3 pointing or sources. 

We examined each object in the H09 list and created a final list of GC candidates from 
the H09 study to include in the combined {HST plus WIYN) radial profile. The final HST 
GC candidate list has 21 objects and includes the "Best" candidates from H09, plus several 
objects that were marked "Bulge" or "Disc" by H09 (to denote that they might be massive 
open clusters in the Bulge or Disk of NGC 891). None of the objects that were noted as 
possible contaminants by H09 (e.g., designated "Star?", "Galaxy?", or "LSB") survived the 
source-by-source review process; in some cases, the objects were clearly galaxies, even in the 
ground-based WIYN images. As part of this same process, we also compared the lists of GC 
candidates from H09 and the WIYN data and investigated the differences between them. 
There were very few overlaps between the WIYN and HST GC candidate lists. Fifteen of the 
21 H09 GC candidates were either in regions of the WIYN images that had been masked out 
(regions in the spiral disk where point sources could not be reliably detected) or located off 
the edge of the WIYN images. Four objects were already identified as WIYN GC candidates. 
Another object had been excluded from the WIYN GC candidate list because it was bright 
(V = 19.3); we decided to add it back to the WIYN list since it was so clearly a GC candidate 
in the HST data. Finally, one object had been excluded from the WIYN GC candidate list 
because it is extended in the WIYN images. 
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The contamination level in the final sample of 21 GC candidates from H09 should be 
minimal; NGC 891 is only ~8 Mpc away so GCs are sometimes resolved in ACS images and 
can usually be distinguished from point-source foreground stars. H09 make use of this and 
shape information to remove stars and galaxies from the sample. The H09 data are very 
deep, reaching 50% completeness at V^29 and /~27— 28, whereas the theoretical GCLF for 
this galaxy should peak at V~22.3 and go to zero at V^~26. Because of this, we assumed 
that the GCLF was fully sampled by the H09 data and that any GCLF coverage correction 
would be zero or negligible. 

5. Results 

5.1. Radial Distributions of the GC Systems 

The final, corrected radial distributions of each galaxy's GC system were produced first 
by assigning the GC candidates from WIYN and the HST studies of G03 or H09 to one of 
a series of concentric annuli, according to their projected radial distances from the galaxy 
center. For the WIYN data, we used annuli of width 1' (as explained in Section I4.2.3p . 
The HST GC candidates were more densely packed and located within a smaller region, 
so were binned more finely, into annuli of width 0.5'. The inner boundary of the first 
annulus was always placed just inside the radial position of the innermost GC candidate. 
An effective area, that excludes missing or masked regions, was calculated for each annulus. 
The number of GCs in each annulus was corrected for contamination and GCLF coverage 
when appropriate, using the corrections discussed in previous sections. We then divided the 
corrected number of GCs by the appropriate effective area to yield the surface density of 
GCs in each annulus. The error on the surface density includes Poisson errors on the number 
of GCs and the number of contaminating objects. The final radial distributions for the GC 
systems of NGC 891 and NGC 4013 are listed in Tables [3] and [Hand plotted in Figures [H] 
and [3 The radii listed in the tables are the average projected radii of the unmasked pixels in 
the corresponding annulus. For each projected radius, the tables list the surface density and 
error, the fraction of the annulus that was observed (excluding missing or masked regions), 
and the source of the measurement (whether it came from the WIYN or HST data). 

We note here that neither of the published HST studies of these galaxies included a 
calculation of the GC system radial profile. G03 stated that low number statistics and 
missing spatial coverage made it too difficult to evaluate the radial distribution of GCs, and 
they estimated the total number of GCs instead by comparing the number of GCs detected 
at various spatial locations around NGC 4013 with the number that would be detected at 
that same location in the Milky Way. They note that this type of analysis implicitly assumes 
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that the spatial distribution of GCs in the target galaxy is similar to that of the Milky Way. 
H09 comment on the "severely limited area coverage" of the ACS study and similarly choose 
to analyze the spatial distribution and population of NGC SQl's GC system by comparing 
the number of GC candidates in their fields with the number that would fall within the same 
projected area around the Milky Way. By combining the data from HST that covers the inner 
few arc minutes of the galaxy with the ground-based WIYN images that extend to larger 
radii, we are able to quantify the radial profile of the GC systems and empirically determine 
how far the systems extend from the host galaxies. This leads to better-determined values 
of the total number and specific frequency of GCs for each system (see the next section). 

The radial distribution of the NGC 4013 GC system (Figure [7]) has the well-behaved, 
regular appearance typical of the GC systems of most of the other spiral galaxies in our 
survey (see RZ03 and R07). The surface density of GCs is at its highest value near the 
galaxy center and then decreases monotonically until it eventually reaches zero within the 
errors, in this case in the r=2.8' (~12 kpc) annulus. Conversely, the profile of NGC 891's GC 
system is not well-behaved. The combined HST plus WIYN profile (Figure [6]) shows positive 
surface density in the first two annuli (from the HST data) inside 1', but then the surface 
density plunges in the third annulus at 1.25', then fiuctuates between zero and positive 
surface density at r~1.5'— 3.4'. All the annuli beyond r~3.4' (~8 kpc) have zero surface 
density within the errors. Without the H09 GC candidates to add to the WIYN sample, we 
would likely have concluded that NGC SQl's GC system was not definitively detected in the 
WIYN images. (The BVR color-color plot in Figure H] and the color-magnitude diagram in 
Figure likewise indicate, at best, a weak WIYN detection of the GC system: there is httle 
or no overdensity of objects within the GC selection box in the color-color diagram, and the 
objects in the color-magnitude diagram are not tightly clustered in the same region of the 
V-{B — R) plane, as they should be.) Likewise, with only the HST data we would not have 
known how far out to integrate the GC system profile in Section 

The corrected radial distributions were fitted with de Vaucouleurs profiles (of the form 
log (Jgc = aO + al r^l^) and power law profiles (of the form log ctgc = aO -|- al log r). For 
NGC 891, the power law had a slightly smaller value (0.9 vs. 1.1 for the de Vaucouleurs 
law); for NGC 4013, the opposite was true, and the de Vaucouleurs was smaller (0.4 vs. 
0.5 for the power law). The best-fit coefficients for both the de Vaucouleurs law and power 
law, for both galaxies' GC systems, are given in Table [51 The function with the smallest 
is plotted in the bottom panels of Figures |6] and [7] and used to calculate the total number of 
GCs in Section [3751 

One of the quantities we evaluate in the wide-field survey is the radial extent of each 
galaxy's GC system, which we define as the point in the radial profile at which the surface 
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density of GCs becomes consistent with zero and remains zero (within the errors) to the outer 
hmit of the data (R07). We tabulated this measure of radial extent for nine survey galaxies 
in R07 and we can now add two more points to the sample. The radial extent for NGC 891 is 
9±3 kpc and for NGC 4013 is 14±5 kpc; the errors include distance errors and uncertainties 
associated with determining the extent from binned radial profile data. The galaxy stellar 



masse s, derived by combining My (Table [T]) with mass-to-light ratios from IZepf fc Ashman 



fll993h . are log(M/M0)=ll.O and 10.9 for NGC 891 and NGC 4013, respectively. In R07, 
we showed that the more massive galaxies in the survey generally have more extended GC 
systems and we derived a quantitative relationship between the GC system radial extent and 
the host galaxy mass. As in R07, we again fit a line and a second-order polynomial to the 
data. The best-fit coefficients of both the line and the curve are consistent within the errors 
with the corresponding coefficients derived in R07. The data and new best-fit second-order 
curve are plotted in Figure [HI The curve has the form 



y = ((49.9 ± 9.2) x^) - ((1080 ± 209) x) + (5860 ± 1190), (1) 

where x is log(M/M0) and y is the radial extent in kiloparsecs. This relationship can be 
used to determine in advance how much spatial coverage is needed to image the GC system 
of a given galaxy. 



5.2. GCs and the Stellar Stream of NGC 4013 



As mentioned in Section NGC 4013 was found bv lMartinez-Delgado et al. I (120091 ) 

to be surrounded by a low-surface-brightness stellar stream, possibly due to a dwarf galaxy 
merging into the larger galaxy's halo. In the process of calculating the asymptotic contam- 
ination correction for the GC system of NGC 4013, we identified 10 GC candidates that 
coincide spatially with the stellar stream and therefore might have originated in the dwarf 
galaxy that is being tidally disrupted. Without velocity measurements for these objects (and 
perhaps even then), we cannot definitively say whether or not the GC candidates are actu- 
ally moving with the stream material. Since only three of the relevant GC candidates have 
V<22.0 and the other seven have 23.4<l^<24.3, velocity measurements would be difficult 
to obtain for the majority of the 10 candidates, given the limits of today's spectroscopic 
facilities. 

We explored various methods for investigating, based on the WIYN images, whether 
there are indeed "stream" GCs in NGC 4013. One straightforward approach is to see whether 
there is an overdensity of GC candidates in the region of the stream loop on the galaxy's 
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northeast side relative to other areas of the GC system. To do this, we defined a rectangular 
region on the WIYN images of NGC 4013 that isolates the part of the tidal stream that 
included the 10 GC candidates that are spatially coincident with the stream. We also defined 
a "control" region with the same total area and mean radial distance from the galaxy center 
and therefore that samples an analogous part of the galaxy's GC system. There are 13 GC 
candidates located within the stream region and 10 in the control region; assuming Poisson 
statistics, these numbers are equivalent. Therefore we cannot say, based on our imaging 
data and a simple number density argument, whether or not there are GCs associated with 
the tidal stream in NGC 4013. This is perhaps not surprising, since it has taken years of 
work - including theoretical modeling, accurate velocity measurements and detailed stellar 
population studies - to show that Sgr dSph in our own galaxy's halo is host to several GCs 
that were once thought to be part of the Milky Way GC system (e.g., Ibata, Gilmore, & 
Irwin 1994, Layden & Sarajedini 2000, Cohen 2004, Carraro et al. 2007). 



5.3. Total Numbers and Specific Frequencies of GCs 

5.3.1. Total Specific Frequency 

We calculated the total number of GCs {Nqc) in each galaxy's system by integrating 
the function that provides the best fit to the final radial distribution. For NGC 891, we 
integrated the best-fitting power law given in Table [51 The inner radius of the integration 
was taken to be 0.2', the inner boundary of the first annulus of the observed radial profile. 
The outer radius of the integration was 3.7', because beyond that radius the surface density 
of GCs is consistent with zero (within the errors) in the rest of the annuli. We had to make 
assumptions about the behavior of the radial profile near the galaxy center, in regions we 
could not observe. For NGC 891, this was the central 0.2' of the GC system, or the central 
~490 pc for m — M =29.61. We calculated the number of GCs that had been missed in 
this inner region given two possibilities: (1) that the profile inside the unobserved region 
continued to follow the best-fit power law to small r, or (2) that the profile fiattened in 
the inner region (i.e., the GC surface density in the unobserved region equaled the value in 
the innermost observed annulus). Adding the number of GCs in the observed portion of the 
system to the number given these two assumptions for the inner region yielded two estimates 
for Ngc'- 53 for the fiat inner profile and 81 if we extended the power law to small r. The 
mean of the two estimates, Nqc = 70±20, is the final value given in Table [61 

For NGC 4013, the de Vaucouleurs law in Table [SI was the best-fitting function to the 
observed radial profile. We integrated the de Vaucouleurs law from 0.08' (the inner boundary 
of the first annulus in the observed profile) to 3.25' (the outer boundary of the 2.8' annulus. 
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which is where the GC surface density reaches zero within the errors and remains zero for 
the rest of the observed profile). For the region inside 0.08', we assumed a fiat radial profile 
or a de Vaucouleurs law all the way to r=0. The results given these two assumptions were 
very close to each other: Nqc = 138 for the flat inner profile and 144 if the de Vaucouleurs 
law was continued to r=0. Again we have taken the average of the two values for Nqc to be 
the final value for the GC system: 140±20 (Table [6]). 

We also calculated luminosity- and mass-normalized specific frequencies for each galaxy's 
GC system. The total number of GCs normalized by the V"-band luminosity of the host 
galaxy is defined as 



S 



N 



(2) 



( iHarris &: van den Bergh Ill98l[ ) and the number of GCs normalized by the host galaxy stellar 
mass is defined as 



T = 



GC 



Mg/109 Mq 



(3) 



(IZepf fc Ashman 1119931 ). We used My from Table [T] and the galaxy mass (designated Mq 
above) from Section ISTTl to calculate Sn and T; specific frequencies are given in Table |6l 

The errors on Nqc given in Table |6] include contributions from the following sources 
of uncertainty: (1) the change in the calculated GCLF coverage depending on the assumed 
intrinsic GCLF dispersion and how the LF data were binned; (2) Poisson errors on the 
number of GCs and contaminating objects; and (3) uncertainties in the number of GCs in 
the inner, unobserved region of the system. The errors on the specific frequencies Sn and 
T also take into account the error on the total galaxy magnitude; we assumed that the 
uncertainty in the extinction-corrected ga laxy magnitude was three times the error on V"^ 
given in RC3 (Ide Vaucouleurs et al. Ill99ll ). Each of these sources of uncertainty was added 
in quadrature to yield the final total error value given in the table. 

Next we examine how our estimates of the number and specific frequency of GCs for 
these galaxies compare to those from past studies. In their HST WFPC2 study of the 
GC system of NGC 4013, G03 estimate the total number of GCs not by constructing and 
integrating the radial profile, but by comparing the number of GC candidates they detect in 
certain regions around the galaxy to the number of GCs that would be detected at analogous 
locations in the Milky Way. Using this method, they calculate that NGC 4013 has Nqc = 
243±51, which is significantly larger than our estimate of 140±20. They combine this with 
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m-M = 31.35 and = -20.83 to derive Sn = 1.1±0.3 and T = 2.2±0.7. They adopt a 
larger distance modulus than we do for NGC 4013, so their inferred galaxy luminosity and 
mass are larger than ours. Our value for Nqc is ~40% lower than their value, but our S^- 
and T values actually agree with the estimates given in G03 because of the larger distance 
modulus they assume. If we instead assume the distance modulus from G03, m — M = 31.35, 
and go through all of the analysis steps to derive final numbers and specific frequencies for 
NGC 4013 (namely, applying a new V magnitude selection to the GC sample, fitting the 
GCLF for both the WlYN and HST data, and calculating and integrating a new radial 
profile), we find Ngc =160±20, Sn = 0.7±0.2, and T = 1.4±0.3. The end result is that, 
even assuming the larger distance to NGC 4013, we derive a smaller total number of GCs in 
the galaxy which leads to correspondingly smaller Sn and T values. 

The GC system of NGC 891 has been studied by van den Bergh & Harris (1982) and 
much more recently by the H09 HST study that we have combined with our WIYN results, 
van den Bergh & Harris (1982) obtained photographic observations with ~1" seeing and 
counted point sources around NGC 891, looking for an overdensity of point-like objects that 
would signal the presence of a GC population. They found "no significant population" of 
GCs in NGC 891 and calculated a specific frequency consistent with zero {Sn ~0.04±0.08); 
they speculated that the weak bulge of NGC 891 may be responsible for its apparent lack 
of GCs. To calculate Nqc, H09 adopt a method similar to that of G03: because their HST 
ACS data also cover only a small portion of the GC system, they compare the number 
of GC candidates they find with the number of Milky Way GCs that would fall within 
a corresponding projected area. From this comparison they conclude that NGC 891 has 
between 80 and 200 GCs and that the GC system "resembles that of the Milky Way rather 
closely" in terms of total number and spatial distribution of GCs. By combining the H09 GC 
candidate list with WIYN imaging data and constructing a radial profile for the GC system, 
we find that NGC 891 has a fairly small total number of GCs (70±20) and small specific 
frequency {Sn = 0.3±0.1, compared to 0.6±0.1 for the Milky Way). Without the HST 
results from H09 to augment our WIYN data, we might have concluded that NGC 891 had 
no detectable GC system as van den Bergh & Harris (1982) did. This galaxy's GC system 
is so poorly populated, and the few GCs that are present are so close to the galaxy disk, 
that typical ground-based capabilities and techniques are simply not adequate to produce a 
definitive detection. On the other hand, combining our ground-based data with the HST 
data allowed us to draw more definitive conclusions than H09 about how many GCs are 
present in NGC 891. 

Figure [3 plots the number of GCs against galaxy Hubble type (top panel) and against 
the log of the galaxy stellar mass (bottom panel) for the spiral galaxies in our wide-field 
survey, plus the Milky Way and M31. Figure [TO] shows Sn versus and T versus the log 
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of the galaxy stellar mass for the same set of galaxies. In R07, we calculated the mean value 
of Ngc^ and the weighted mean values of Sjq and T, for the spiral galaxies in our survey plus 
the Milky Way and M31. These were Nqc = 170±40, Sn = 0.8±0.2, and T = 1.4±0.3. The 
values we calculate for NGC 4013's GC system are very much in line with these mean values, 
but the Ngc, Sn, and T values for NGC 891 are comparatively low. NGC 891 is frequently 
described in the literature as being similar to the Milky Way because of its morphological 
type, luminosity, and molecular gas distribution (e.g., Sandage 1961, van der Kruit 1984, 
Scoville et al. 1993), but its GC system is about half as populous as the Milky Way's. One 
possibility we considered is whether this galaxy has very few GCs because, despite being 
similar to the Milky Way in appearance, it is actually much less massive than our Galaxy. 
As explained in Section EH we calculate log fM/M^T,) = 1 1.0 for NGC 891. The Milky Way 



has My = —21.3 and is an Sbc spiral galaxy ( iHarris Ill99ll ): assuming the M/Ly values from 



Zepf fc Ashman I fll993h . this yields a mass of log(M/MQ) = 11.2. The CO + HI rotation 
curve of NGC 891 is similar to that of the Milky Way, with a rotation velocity that rises 
steeply near the gal axy nucleus to a maximum velocity of 250 km s~^ at ~3 kpc out to 



15 kpc (jSofue 1 119961 ) ■ suggesting this galaxy and the Milky Way probably have comparable 
mass distributions. In any case, the new mean value for Ngc for spiral galaxies in Figure M 
is 150±30. The weighted mean values for Sn and T for the galaxies in Figure ITO] are 0.5±0.1 
and 1.1±0.2, respectively. These updated values agree within the errors with results from 
Chandar et al. (2004), who studied five spiral galaxies and found a mean Sn of 0.5±0.1 and 
T of 1.3 ±0.2. However, our new mean values are slightly smaller than the average values 
from the G03 study of five Sab— Sc spiral galaxies; G03 found a mean Sjy of 0.96±0.26 and 
Tof 2.0±0.5. 

The sample of spirals we have observed with the wide-field survey is still fairly small, not 
all morphological types are sampled, and we have more Sb spiral galaxies in the survey than 
any other type. Therefore we do not wish to over-interpret the significance of the results. 
However, based on this small sample, there does seem to be modest variation in the number 
of GCs and in the GC specific frequency for galaxies with similar properties. For example, 
the Sb galaxy NGC 4157 and the Sab galaxy NGC 7814 have almost the same stellar mass 
(log(M/M0 = 10.9), but the number of GCs (80±20 and 170±30, respectively) differs by 
a factor of two. NGC 891, NGC 2683, NGC 4157, and NGC 4013 are all classified as Sb 
spiral galaxies, but the number of GCs they host ranges from 70 to ^200. The bottom line 
seems to be that in the overall sample, there may be a weak trend of increasing Nqc with 
increasing galaxy mass (bottom panel of Figure [9]) , but there is enough variation in Nqc at 
each galaxy mass, magnitude, and/or Hubble type to produce scatter in the T versus galaxy 
stellar mass and Sj^ versus galaxy magnitude relationships seen in Figure [TDl 
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5.3.2. Specific Frequency of Blue (Metal-Poor) GCs 



Many giant galaxies - including the two with arguably the most thoroughly studied GC 
systems, the Milky Way and M31 - have been shown to possess GC systems with at lea st two 



populations: a blue, metal-poor population and a redder, more metal-rich populations (IZinn 



19851 : IZepf fc Ashman 



20021: 



Kundu fc Zepf 



1993 



2001 



Barmbv et al. II2OOOI: iKundu fc Whitmore II2OOII : IPerrett et al. 



3; Istrader et alTlboOTh . Various galax y formation models h ave 



either predicted the presence of these multiple populations of GCs (lAshman fc Zepf Ill992) or 



sought to explain them in the context o f the galaxies' assembly histories (IForbes et al. 



1997 



Cote et al. Ill998l : iBeasley et al. II2OO2I ). The metal-poor population of GCs is of particular 
interest because its low metallicity implies that it represents the first generation of GCs 
created in the early stages of the host galaxies' formation histories. 

In previous papers from the wide- field survey (R05, ROT), we estimated the number of 
blue (metal-poor) GCs normalized by the galaxy mass, or Tbiue^ 



^'"^ - Mg/109 Mq ^ > 

Here, Ai'G'c(blue) is the number of blue GCs and Mq is the stellar mass of the host galaxy, 
calculated as described in Section 15.11 In R05 and R07, we found that Tbiuo for the spiral 
galaxies in the survey is significantly smaller than Tbiue for the massive cluster ellipticals. This 
leads to the conclusion that simply merging the GC systems of two or more late-type galaxies 
as we see them today cannot account for the comparatively large blue GC populations in 
some giant ellipticals. We also found a rough trend of increasing Tbiue values with increasing 
galaxy mass. 

We calculated Tbiue for NGC 891 and NGC 4013 to see how they compare to the values for 
the other survey galaxies, and whether they fit within this overall trend. We first construct 
a sample of GC candidates with magnitude completeness of at least 90% in all the filters; we 
refer to this sample of candidates as the "complete sample" . We then determine the fraction 
of GC candidates in the complete sample that are bluer than some pre-determined color. 
The color criterion we use is based on where the usual separation between the blue and red 
GC populations is located in elliptical galaxies; m. B — R this occurs at ~1.23 (RZOl, RZ04) 
and in — / it occurs at ~1.0 (Kundu & Whitmore 2001). 

For NGC 891, we used the WIYN data to construct a complete sample of GC candidates 
and found that seven of the 26 candidates (27%) in the complete sample had B — R bluer 
than 1.23. Because we had so few objects in the complete sample from WIYN, we also 
decided to construct a complete sample of GC candidates using the HST list from H09. 
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The complete sample from the H09 data included 21 GC candidates; 13 of these (63%) had 
V — I < 1.0. We multiplied each of these blue fractions times the T value for the galaxy's 
GC system to calculate Tbiue and then took the mean of the two Tbiue values as the final 
value. We also assigned a generous error that incorporates a wide range of possible values 
for Tbiue, given that working with very small samples of GC candidates makes the entire 
process of the Tbiue calculation for these spiral galaxies very uncertain. The final Tbiue value 
for NGC 891 is 0.3±0.2. For NGC 4013, the complete sample from the WIYN data included 
51 GC candidates; of these, 38 (75%) have B — R < 1.23. Multiplying this fraction times the 
calculated T value for NGC 4013's GC system yields a Tbiue value of 1.4±0.3. These values 
are listed, along with the other total numbers and specific frequencies, in Table [61 

In Figure [TTl we present T^ue values for 16 galaxies, plotted versus the galaxies' stellar 
masses. Note that this is the same Tbiue figure presented in R07, but with NGC 891 and 
NGC 4013 added. The figure now includes: 11 spiral, SO, and elliptical galaxies from our 
wide- field survey (RZOl, RZ03, RZ04, R05, R07); three early-type galaxies from other multi- 
color, wide-field CCD studies in the literature (Forbes et al. 2001, Gomez & Richtler 2004, 
Harris et al. 2004); and the Milky Way and M31 (Zinn 1985, Ashman & Zepf 1998, Barmby 
et al. 2000, Perrett et al. 2002). Circles designate giant elliptical galaxies in clusters, squares 
are used for early-type galaxies in the field, and triangles designate spiral galaxies in the field 
or in groups. Tbiue for the Milky Way and M31 are the triangles with relatively small error 
bars: T^ue = 0.9±0.2 and 1.2±0.3, respectively. The T^ue value for NGC 4013 (1.4±0.3) 
is very close to the mean derived for the other spiral galaxies without the two new points 
(1.2±0.1). Tbiuc for NGC 891, however, is relatively low, at only 0.3±0.2; this is comparable 
to the value we found for the spiral galaxy NGC 7331, 0.4±0.3 (R07). 

The general trend identified in R05 — i.e., larger numbers of blue, metal-poor GCs 
per galaxy stellar mass for more massive galaxies — is still present in this latest version of 
the Tbiue figure, although (as before) there is a fair amount of scatter in the relationship. 
In R05 we suggested that the general trend appears consistent with a scenario in which 
galaxies and their GC systems are formed in a biased, hierarchical assembly process. As 
explained in Santos (2003), the idea is that the formation of giant galaxies began at high 
redshift {z > 20) as gaseous protogalactic building blocks collided and merged together to 
form larger galaxies. Merging of gaseous protogalaxies triggers the formation of GCs (in a 
similar way as envisioned by Ashman & Zepf 1992). As galaxies are built up hierarchically 
over time, they increase both their total stellar mass and their GC population. Massive 
galaxies that are located in parts of the universe that have higher density (e.g., areas of 
the universe where galaxy clusters will eventually be located, at 2; ~ 0) begin this assembly 
process earliest, before less-massive galaxies that are located in lower-density environs. In the 
Santos (2003) picture, galaxy assembly and GC formation are temporarily suppressed when 
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reionization occurs in the universe. The energy imparted to the intergalactic medium slows 
down or stops protogalactic building blocks from merging and forming stars, perhaps for a 
period of 1 Gyr. Stars continue to evolve during this time, so that when the assembly and 
associated GC formation process resumes, the gas has been enriched and any GCs formed 
after the suppression period are metal-rich by comparison. The result of such a scenario 
is that today's more massive galaxies should have higher numbers of metal-poor GCs per 
galaxy mass than less massive galaxies, simply because they started the assembly process 
first; a larger percentage of their baryonic mass was in place to participate in the formation of 
the first generation of GCs before the temporary suppression occurred. This type of scenario 
also provides a natural explanation for the gap in metallicity between GC subpopulations in 
the Milky Way and other giant galaxies. 

The curves in Fi gure [TT] were made by G. Bryan (private communication) using a Press- 



Schecter calculation ( jPress fc Schechter Ill974l ) first described in R07. The curves show how 



the slope of the Tbiue trend would be expected to change as the redshift of suppression (i.e., 
the redshift at which galaxy assembly and GC formation are temporarily stopped) changes, 
for three choices of redshift: 2; = 7, 11, or 15. Bryan's calculation assumes that GCs form 
within gaseous building blocks with masses of at least 10^ solar masses and that the number 
of GCs in a given galaxy is proportional to the fraction of the galaxy mass that has assembled 
by the suppression redshift. The current data do not seem to match the z = 7 curve, but 
instead suggest a larger suppression redshift of z~ll — 15. (Given that Santos (2003) suggests 
that reionization is the mechanism that halts the formation of the first generation of GCs, 
it is interesting to note that the z~ll— 15 curve matches well with recent estimates for the 
reionization redshift from WMAP (Alvarez et al. 2006)). 

The biased, hierarchical scenario outlined here is not the only way that the trend in Tbiue 
seen in Figure [TT] could arise; other possible explanations (e.g., differences in the efficiency 
of dynamical destruction of GCs in different types of galaxies, a trend in the relationship 
between galaxy luminosity and galaxy stellar mass that would affect the log(M/M0) values 
in the figure) are discussed in R05 and R07. Our intent here is simply to show how the 
Tbiue and stellar mass values of NGC 891 and NGC 4013 compare with the values for the 
other survey galaxies. The figure — and our understanding of the relationship between Tbiue 
and galaxy and GC system formation and evolution — would certainly benefit from having 
additional well-determined Tbiue values at the high-mass end (where we only have values for 
two Virgo cluster ellipticals) and in the mo derate- luminosity region (log(M/M0) ^11.5— 12). 
Accordingly, we have an ongoing effort to observe more galaxies with multi-color wide-field 
imaging and fill in the gap in the relevant mass range. 
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6. Summary and Conclusions 

We have combined BVR imaging observations from the 9.6'x9.6' Minimosaic imager on 
the WIYN 3.5-m telescope with archival and published data from HST to investigate the 
GC system properties of two spiral galaxies, NGC 891 and NGC 4013. The main findings of 
this work are as follows: 



1. We used both the WIYN data and the HST data to construct radial distributions of the 
GC systems of the galaxies, corrected to account for missing areal coverage, the portion 
of the GCLF we could not observe given our detection limits, and contamination from 
foreground stars and background galaxies. NGC 891's GC system radial profile extends 
to 9±3 kpc and NGC 4013's extends to 14±5 kpc from the centers of the host galaxies; 
beyond those projected radii, the surface density of GCs is consistent with zero within 
the errors. 

2. We fitted the radial distributions with de Vaucouleurs and power laws and integrated 
the best-fitting function to derive the total number of GCs. We find Nqc = 70±20 
for NGC 891 and 140±20 for NGC 4013. Despite the fact that NGC 891 is often 
considered a Milky Way analog, we estimate that it has less than half as many GCs 
as the Galaxy. Its specific frequency is 0.3±0.1 and its mass-normalized number of 
GCs, T, is 0.6±0.3, com pared to = 0.6±0.1 and T = 1.3±0.2 for the Milky Way 



flAshman k Zepf Ill998f ). We find that NGC 4013 has Ngc = 140±20, Sn = 1.0±0.2, 



and T = 1.9±0.5. 

3. We combine the results from this paper with results from our ongoing wide-field GC 
system survey and calculate the average total number and specific frequency of GCs 
for spiral galaxies in the survey plus the Milky Way and M31. We find Nqc = 150±30, 
S'tv = 0.5±0.1, and T = 1.1±0.2. These values are consistent with values from our 
most recent survey paper (R07) and from an analysis of five spiral galaxy GC systems 
by Chandar et al. (2004) but slightly smaller than those from a study of several spiral 
galaxies by Goudfrooij et al. (2003). Plots of number and specific frequency of GCs 
versus galaxy Hubble type, absolute magnitude, and stellar mass for our sample show 
no strong trends and a modest amount of scatter in GC system properties for otherwise 
similar host galaxies. 

4. We derive the number of blue, metal-poor GCs normalized by the galaxy stellar mass 
(^biue) for both NGC 891 and NGC 4013 and find them generally consistent with a 
trend of increasing specific frequency of blue GCs with increasing galaxy mass. As 
was discussed in previous survey papers (e.g., R05, R07), the Tbiue trend may be a 
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consequence of biased, hierarchical galaxy formation. In such a scenario, a larger 
fraction of the mass of high-mass galaxies is in place when the first generation of GCs 
is formed; the consequence is that higher-mass giant galaxies end up with larger Ty^i^e 
values than their lower-mass counterparts. 

Finally, a key conclusion of this work is that, although this is not always possible, the 
most effective way to measure the global properties of the entire system of GCs around 
giant galaxies in the local universe is to combine good-quality wide-field imaging with HST 
observations. The ground-based images allow one to trace a GC system to its outer radial 
extent, while the HST data help to characterize contamination and probe further into the 
galaxy disk. One can then use both data sets to carefully construct a list of GC candidates 
and a radial profile for the system and derive more accurate global quantities that describe 
the GC population. 
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Fig. 1. — Digitized Sky Survey image of NGC 891 with the location of the WIYN point- 
ing (large box) superimposed. The three smaller boxes marked HI, H2, and H3 show the 
locations of the HST ACS fields analyzed by Harris et al. (2009) and discussed in Section 1131 
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Fig. 2. — Digitized Sky Survey image of NGC 4013 with the location of the WIYN pointing 
(large box) superimposed. The HST WFPC2 pointing analyzed by Goudfrooij et al. (2003) 
and discussed in Section 14.41 is also marked. 
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0.5 1 1.5 2 

(B-V)o 

Fig. 3. — Color selection of GC candidates in NGC 4013. Open squares are the 488 point 
sources detected in all three filters; filled circles are the final set of 83 GC candidates. Because 
V magnitude criteria were imposed, not all objects within the color selection box are GC 
candidates. A reddening vector of length Ay = 1 mag appears in the upper left-hand corner. 
The box shows the boundaries of the color selection criteria. The curves show where galaxies 
with different morphological types, located at redshifts from 2; = to 0.7, would lie in this 
color-color plane. 
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Fig. 4. — Color selection of GC candidates in NGC 891. Open squares are the 490 point 
sources detected in all three filters; filled circles are the final set of 43 GC candidates. 
Because V magnitude criteria were imposed, not all objects within the color selection box 
are GC candidates. The other features of the plot (reddening vector, galaxy curves) are as 
in Figure [3l 




Fig. 5. — V magnitude versus B — R color for the final set of GC candidates in NGC 891 
and NGC 4013. The magnitudes and colors have been corrected for Galactic extinction in 
the direction of the galaxy fields, but not for absorption internal to the galaxies themselves. 
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Fig. 6. — Radial distribution of GCs in NGC 891, plotted as surface density vs. projected 
radial distance, r (top), and as the log of the surface density vs. log of r (bottom). The 
horizontal dotted line in the top panel marks zero surface density. The dashed line in the 
bottom panel is the best-fit power law. The data have been corrected for contamination, 
areal coverage, and magnitude incompleteness, as described in Section 15.11 
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Fig. 7. — Radial distribution of GCs in NGC 4013, plotted in the same way as in Figure [61 
The dashed line in the bottom panel is the best-fit de Vaucouleurs law. 
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Fig. 8. — Radial extent of the GC system in kiloparsecs versus the log of the galaxy stellar 
mass in solar masses for 11 elliptical, SO, and spiral galaxies included in our wide-field GG 
system survey to date. Open triangles are the points from R07; filled circles are two new 
points from the current paper. The best-fit second-order polynomial (which is the same, 
within the errors, as the one given in R07) is shown as a dashed line. 
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Fig. 9. — Number of GCs plotted vs. Hubble type (top) and galaxy stellar mass (bottom). 
Values for the spiral galaxies analyzed for our wide-field survey are plotted with filled circles 
(the two galaxies from the current paper) or open circles (four galaxies from R07 and one 
from RZ03). Values for the Milky Way (Ashman & Zepf 1998) and M31 (Ashman & Zepf 
1998, Barmby et al. 2000; larger error bars) are plotted with open stars. 
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Fig. 10. — vs. y-band galaxy luminosity (top) and T vs. galaxy stellar mass (bottom). 
Symbols are the same as in Figure O 
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Fig. 11. — Galaxy- mass-normalized specific frequency of blue (metal-poor) GCs plotted ver- 
sus the log of the galaxy stellar mass for sixteen galaxies from our widc-ficid survey and 
from other multi-color wide-field CCD imaging studies in the literature. Circles represent 
early-type galaxies in clusters; squares indicate early-type field galaxies; and triangles indi- 
cate field spiral galaxies. Open symbols denote the two new values from the current study. 
The curves (courtesy G. Bryan) show the slope of the Tbiue vs. galaxy stellar mass relation 
depending on whether metal-poor GCs were formed at 2; > 7, 11, or 15 (see the text for 
details) . 
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Table 1. Basic Properties of the Target Galaxies 



Name 


Type 


Inclination 


VHelio 


m- M 


Distance 








(deg) 


(km s~^) 




(Mpc) 




NGC 891 


Sb 


84 


831±1 


29.61 


8.36 


-20.79 


NGC 4013 


Sb 


90 


528±4 


30.90 


15.1 


-20.38 



Note. — Morphological types are from the Hubble Atlas of Galaxies (Sandage 
1961) for NGC 891 and from RCl (de Vaucouleurs & dc Vaucouleurs 1964) for 
NGC 4013. Inclination and heliocentric radial velocity arc from Vcrheijen & 
Sancisi (2001) for NGC 4013. Heliocentric radial velocity for NGC 891 is from 
RC3 (de Vaucouleurs et al. 1991) and inclination is from TuUy (1988). Distance to 
NGC 891 is from Tonry et al. (2001) (surface brightness fluctuations). Distance 
to NGC 4013 is from combining the recession velocity with respect to the 3 K 
cosmic microwave background from the NASA Extragalactic Database with Hq 
= 70 km s~^ Mpc~^. Total absolute magnitudes are from combining Vj? from 
RC3 with m-M. 



Table 2. WIYN Minimosaic Observations of the Target 

Galaxies 



Galaxy 


Date 


Exposure Times ( 








B V 


R 


NGC 891 


2001 Jan 


4 X 2100 3 X 2000 3 


X 1500 


NGC 4013 


2001 Jan 


3 X 2100 1 X 2000 3 


X 1500 




2009 Mar 


4 X 2000 





-42 - 



Table 3. Corrected Radial Profile of the GC System of NGC 891 



Radius Surface Density Fractional Coverage Source 
(arcmin) (arcmin~^) 



n 48 


9 62 + 3 21 


66 


HST 


u.yt) 


4. DC) ± i.oy 




WQT' 

rlo i 


1.25 


0.64 ± 0.57 


0.71 


WIYN 


1.46 


0.61 ± 0.61 


0.36 


HST 


1.97 


0.43 ± 0.43 


0.38 


HST 


2.17 


0.91 ± 0.45 


0.77 


WIYN 


2.46 


0.00 ± 0.00 


0.42 


HST 


2.95 


0.52 ± 0.37 


0.41 


HST 


3.16 


0.43 ± 0.33 


0.64 


WIYN 


3.44 


0.37 ± 0.37 


0.25 


HST 


3.94 


0.00 ± 0.00 


0.16 


HST 


4.15 


0.10 ± 0.28 


0.47 


WIYN 


4.45 


0.57 ± 0.57 


0.13 


HST 


5.16 


0.25 ± 0.29 


0.43 


WIYN 


6.12 


-0.01 ± 0.25 


0.32 


WIYN 


7.12 


0.09 ± 0.30 


0.22 


WIYN 


8.04 


-0.09 ± 0.40 


0.07 


WIYN 



Note. — For completeness, we list all points in the radial profile 
made using the WIYN + HST GC candidate samples; when no GC 
candidates appear within a particular annulus, the surface density 
is written as "O.OOiO.OO". Negative surface densities sometimes 
occur because a contamination correction has been applied to each 
annulus. 
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Table 4. Corrected Radial Profile of the GC System of 

NGC 4013 



Radius 


Surface Density 


Fractional Coverage 


Source 


(arcmin) 


(arcmin"^) 






0.39 


37.30 ± 6.93 


1.00 


HST 


0.82 


9.90 ± 2.65 


0.72 


HST 


0.90 


9.70 ± 2.31 


0.72 


WlYN 


1.34 


7.07 ± 2.67 


0.32 


HST 


1.81 


2.22 ± 1.08 


0.89 


WIYN 


2.78 


0.90 ± 0.93 


0.96 


WlYN 


3.77 


0.37 ± 0.89 


0.96 


WIYN 


4.68 


0.36 ± 0.90 


0.66 


WIYN 



Table 5. Coefficients from Fitting Radial Proffie Data 



de Vaucouleurs Law Power Law 



Galaxy aO al aO al 



NGC 891 3.02 ± 0.33 -2.53 ± 0.29 0.45 ± 0.07 -1.62 ± 0.18 
NGC 4013 4.02 ± 0.30 -3.11 ± 0.32 0.88 ± 0.05 -1.73 ± 0.18 



Table 6. Total Numbers and Specific Frequencies of 
the GC Systems 



Galaxy 


Ngc 




T 


^blue 


NGC 891 


70±20 


0.3±0.1 


0.6±0.3 


0.3±0.2 


NGC 4013 


140±20 


1.0±0.2 


1.9±0.5 


1.4±0.3 



